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1.0  Introduction 

A  Specialists  Meeting  on  Atomization  and  Non-dllute  Spray  Behavior  was  held 
on  March  20  and  21,  1985  at  SandU  National  Laboratories.  The  meeting  was 
jointly  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  the  Army 
Research  Office  and  Sandia  National  Laboratories.  The  convening  of  the 
meeting  was  motivated  by  the  perception  that  the  areas  of  atomization  and 
non-dilute  spray  behavior  were  technologically  Important,  poorly  understood 
and  little  studied.  It  was  the  expectation  of  the  sponsors  that,  by  providing 
a  forum  for  experts  to  review  the  current  state  of  knowledge  and  identify 
critical  deficiencies,  it  would  be  possible  to  establish  directions  and  priori 
ties  for  increasing  the  research  in  these  areas. 

Attendance  at  the  meeting  was  by  Invitation  only.  This  was  done  to  keep  the 
number  of  attendees  snail  enough  to  promote  active  discussion  and  participa¬ 
tion  from  all.  An  attempt  was  made  to  include  researchers  who  were  representa 
tive  of  current  activity  in  spray  research,  as  well  as  those  working  in  poten¬ 
tially  relevant  areas.  The  List  of  Attendees  Is  Included  as  Appendix  A. 

The  meeting  was  organized  into  six  sessions.  The  first  four.  Dilute  Spray 
Behavior,  Non-dilute  Spray  Behavior,  Atomization,  and  Measurement  Techniques, 
were  designed  to  review  past  and  current  research  and  provide  a  framework  for 
the  last  two  sessions.  These  were  a  General  Discussion  of  the  material  pre¬ 
sented  in  the  four  reviews  and  final  session  devoted  to  developing  a  Summary 
and  Research  Recommendations  resulting  from  the  meeting.  In  the  subsequent 
sections  of  this  report,  each  of  the  review  sessions  is  summarized,  followed 
by  a  synopsis  of  the  recommended  research  needs  in  the  field.  The  Figures 
used  by  the  four  invited  speakers  are  Included  as  Appendices  B  through  E. 

These  Figures  ire  referenced  In  the  following  summary  sections. 

2.0  Dilute  Sprays 

The  summary  of  Dilute  Spray  Behavior  was  given  by  Professor  Gerard  Faeth 
of  Pennsylvania  State  University.  The  lecture  was  Intended  to  preface  the 
areas  of  major  concern  for  the  meeting  by  reviewing  results  In  a  spray  regime 
in  which  research  directions  were  relatively  well  defined  and  not  In  need  of 
further  stimulation.  In  fact,  even  this  presentation  caused  considerable 
discussion  among  the  participants. 

The  fourth  figure  (see  Appendix  B)  of  Professor  Faeth's  presentation  defines 
dilute  sprays  based  on  behavioral  limitations.  This  figure  provoked  a  strong 
controversial  response.  It  became  clear  that  the  boundary  between  dilute  and 
non-dilute  behavior  itself  Is  a  significant  unknown  and  may  be  coupled  to  the 
spray  phenomenon  of  Interest.  Thus,  the  boundary  for  droplet  interactions 
Influencing  evaporation  or  group  mode  combustion  may  differ  from  that  for 
turbulence  modulation.  Figure  4  then  may  be  regarded  as  a  general  rule  of 
thumb,  reflecting  conditions  which  may  be  neither  necessary  nor  sufficient  for 
dilute  spray  behavior. 

Figures  8-11  Identify  three  approaches  to  modeling  spray  behavior.  These 
approaches  reflect  contemporary  capability  for  predicting  turbulent  flow 
behavior  In  which  drop  sizes  are  much  smaller  than  grid  resolution  of  compu¬ 
ters.  Therefore,  spray  behavior  has  been  modeled  as  a  subgrid  phenomenon 
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using  Reynolds  or  Favre-averaged  differential  transport  equations  with  various 
levels  of  closure.  The  results  given  In  Professor  Faeth's  presentation  all 
were  calculated  with  the  k*eps11on  closure  model. 

The  first  of  the  three  models  for  particle-laden  flows,  illustrated  in  the 
eighth  presentation  figure,  is  the  locally  homogeneous  flow  (LHF)  model.  In 
this  model  discrete  particle  phenomena  are  ignored,  and  the  particulate  is 
modeled  as  a  jet  of  a  second  fluid  continuum  injected  into  the  ambient  fluid. 
The  model  is  rigorously  vali'^  only  for  infinitely  small  particles  at  low 
volume  fractions. 

The  second  model  discussed  by  Professor  Faeth  and  illustrated  in  Figure  9  is 
the  discrete  separated  flow  (OSF)  model.  In  this  model  a  particle  laden  flow 
is  sfmul'ated  by  a  compu'ta'tTona  1 1  y  tractable  number  of  discrete  particles. 

Mass,  momentum  and  energy  exchanges  between  the  particle  and  the  fluid  ambient 
are  influenced  by  the  mean  velocity  of  the  fluid  flow  and  not  by  ambient 
turbulence  behavior. 

The  last  of  the  three  models,  shown  in  the  tenth  figure,  is  the  stochastic 
separated  flow  (SSF)  model.  It  differs  from  the  OSF  class  of  models  in  that 
amb i ent  tu rbuTence  does  influence  particle  behavior  through  a  random  sampling 
of  turbulent  eddy  behavior. 

Figures  14-19,  20-22,  23-25,  26-28  and  29-39  show,  respectively,  comparisons 
between  the  predictions  of  the  three  models  and  experimental  measurements  for 
the  particle-laden  flows  listed  in  Figure  13.  Note  that  solid  particles  in 
air  and  gaseous  bubbles  in  liquid  were  tested  along  with  liquid  sprays  in  air. 
The  SSF  model  predictions  generally  are  quite  accurate  and  outperform  the 
other  models.  Some  inaccuracies  are  noteworthy  for  SSF  predictions  of  radial 
dispersion,  but  these  errors  are  attributed  to  the  inability  of  the  k-epsilon 
closure  model  to  account  for  anisotropic  turbulence.  These  predictions 
reflect  contemporary  spray  prediction  capability. 

Outstanding  research  issues  for  dilute  spray  behavior  are  summarized  in  Figure 
57.  Figures  42-56  reflect  the  current  state  of  understanding  for  these 
Issues,  including  research  programs  which  are  currently  active.  These  Figures 
suggest  that  significant  research  remains  for  single  droplet  and  dilute  spray 
behavior.  A  strong  coupling  exists  to  the  atomization  and  nondilute  phenome¬ 
na,  which  Initialize  the  aerothermochemical  behavior  of  dilute  sprays. 

3.0  Non-dilute  Sprays 

The  review  of  non-dilute  sprays  was  given  by  Professor  C.  K.  Law  of  the  Univer¬ 
sity  of  California,  Davis.  Suggested  references  were  the  review  article  by 
Faeth  (Reference  1)  and  the  paper  by  O'Rourke  and  Bracco  (Reference  6). 

As  shown  in  Figure  3  (see  Appendix  C),  the  non-dilute  or  dense  spray  is  a  tran¬ 
sition  state  between  the  atomization  region  and  the  dilute  spray.  In  this 
region  the  droplet  density  is  high  so  that  droplet-droplet  Interactions  are 
Important.  Further,  the  spray  is  not  fully  mixed  with  the  ambient.  Thus 
there  are  substantial  radial  temperature  and  velocity  gradients  and  droplets 
are  moving  at  moderate  Reynolds  numbers.  The  extent  of  the  non-dilute  spray 
region  is  primarily  governed  by  the  rate  of  ambient  entrainment  and  mixing. 

The  spray  is  dilute  when  the  inter-droplet  spacing  has  increased  to  the  point 
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that  droplets  may  be  considered  as  interacting  with  the  intervening  fluid 
instead  of  each  other.  A  code  developed  for  the  non-dilute  spray  region 
should  be  capable  of  treating  the  evolution  of  the  initial  droplet  size  and 
momentum  distributions,  as  given  by  the  appropriate  atomization  model.  At 
present,  this  initial  input  is  largely  empirical. 

Non-dilute  sprays  have  been  regarded  as  presenting  a  more  difficult  challenge 
for  analysis  than  dilute  sprays  because  of  phenomena  described  in  the  current 
summary.  Dr.  Law  suggested  that  opportunities  for  simplifications  in  analyses 
in  non-dilute  spra>s  also  should  be  explored  (see  Figures  6  and  7).  For 
example,  the  interior  of  the  spray  may  be  relatively  cool  and  saturated  with 
vapor  so  that  mass  transfer  and  change  of  phase  effects  can  be  ignored. 

While  there  may  be  differences  in  approach  to  the  detailed  modeling,  there 
appears  to  be  agreement  on  the  important  processes  which  must  be  treated  in 
a  description  rf  non-dilute  sprays.  These  are  droplet  collision,  droplet 
transport,  and  two-phase  flow  interactions. 

Droplet  collisions  may  result  in  coalescence  (see  Figures  14-20),  with  one 
resultant  large  droplet,  partial  coalescence  followed  by  separation,  resulting 
in  the  change  in  sizes  of  collision  partners  and/or  the  creation  of  new  drop¬ 
lets,  or  inelastic  recoil,  with  no  change  in  size.  Coalescence  is  governed  by 
the  rotational  energy  of  the  collision  complex.  If  the  roLational  energy  of 
the  complex  exceeds  the  increased  surface  energy  necessary  to  form  two  drop¬ 
lets,  the  complex  will  be  unstable  and  separate.  There  is  little  quantitative 
information  on  the  other  possible  processes.  Therefore,  the  inclusion  of 
droplet-droplet  interactions  involves  the  development  of  approximate,  largely 
intuitive  models  for  collision  effects. 

Droplet  transport  processes  in  non-dilute  sprays  must  be  modeled  over  the 
range  spanning  fluidized  bed  to  non-interacting  conditions.  In  addressing 
these  problems,  Dr.  Law  presented  an  extensive  review  of  vaporization  pro¬ 
cesses  relevant  to  non-dilute  sprays.  This  discussion  is  summarized  in 
Figures  9-15.  Approximate  relations  have  been  developed  for  drag  coeffici¬ 
ents,  heat  and  mass  transport  rates  in  dense  flows.  However,  these  are 
largely  based  on  data  obtained  for  particulate  or  low  vapor  pressure  fluids. 
Further  research  is  needed  to  develop  appropriate  relations  for  vaporizing 
droplets.  Dense  spray  mass  transport  becomes  particularly  important  for 
combusting  flows.  The  dynamics  of  group,  as  opposed  to  individual  droplet, 
combustion  are  governed  b;'  the  fuel  vaporization  and  transport  rate  from  the 
fuel  rich,  dense  spray  to  the  outer,  oxygen-rich  mixing  zone. 

In  the  non-dilute  spray  region,  the  ambient  gas  is  being  accelerated  and 
entrained  by  the  high  velocity  droplet  stream.  Turbulence  production  is 
reduced  due  to  particle  inertia  which  also  modifies  turbulent  transport  (see 
Figure  21). 

4.0  Atomization 

The  review  of  jet  atomization  was  presented  by  Dr.  Rolf  Reitz  of  General 
Motors  Research  Laboratory  and  followed  the  general  outline  of  Reference  11. 
While  there  appears  to  be  no  precise  definition  of  atomization,  for  the  pur¬ 
poses  of  this  discussion,  the  term  will  be  associated  with  that  region  of  the 
injected  fuel  spray  in  which  the  fuel  is  reduced  to  droplets  of  size  much  less 


than  the  nozzle  orifice.  Thus,  the  atomization  region  can  be  considered  as 
extending  from  an  undefined  point  upstream  of  the  nozzle  orifice  to  the  non- 
dllute  region  of  the  spray. 

With  other  conditions  held  fixed,  as  shown  In  Figure  3  (see  Appendix  0),  It 
has  been  possible  to  Identify  four  distinct  regimes  of  jet  breakup  which  are 
a  function  of  jet  velocity.  At  very  low  speeds,  the  jet  Is  unstable  with 
respect  to  small  perturbations  and  forms  droplets  of  a  size  equal  to  or  larger 
than  the  jet  (Rayllegh  Instability  regime).  At  slightly  higher  velocities, 
the  jet  experiences  a  low  shear  interaction  with  the  amolent  gas  and  jet 
breakup  Is  augmented  by  the  aerodynamic  Interaction  between  the  liquid  and  gas 
(First  wind-induced  breakup  regime).  The  droplet  size  Is  still  on  the  order 
of  the  jet  diameter.  At  higher  velocities,  jet  breakup  Is  due  to  the  growth 
of  unstable  surface  waves  produced  by  the  relative  motion  between  the  jet  and 
ambient  gas  (Second  wind-induced  breakup  regime).  Since  the  wavelength  of 
these  disturbances  Is  short,  the  resulting  droplets  are  much  less  than  the  jet 
diameter.  With  the  possible  exception  of  the  second  wind-induced  breakup 
regime,  the  above  processes  may  have  little  to  do  with  the  formation  of  spray 
from  practical  fuel  nozzles,  since  all  three  Involve  breakup  at  a  considerable 
distance  downstream  of  the  nozzle  exit.  In  contrast,  both  single  and  dual 
fluid  (air  assist)  atomizers,  when  operating  at  their  design  points,  appear  to 
atomize  the  fluid  almost  Instantaneously  upor  emergence  of  the  jet  from  the 
nozzle.  This  final  regime  of  jet  breakup  has  been  termed  “atomization'  and  is 
characterized  by  high  jet  velocities  and/or  high  relative  gas-liquid  veloci¬ 
ties. 

Also  in  contrast  to  the  other  three  regimes,  there  Is  a  lack  of  consensus  on 
the  mechanisms  operating  In  the  atomization  regime.  Due  to  small  dimensions, 
the  high  liquid  densities  and  other  impediments  to  experimental  Investigation, 
there  Is  little  direct  knowledge  about  the  stat>  of  the  liquid  and  gas  In  this 
region.  Experimental  studies  have  concentrated  on  elucidating  the  variables 
controlling  the  spray  shape  and  droplet  size  distribution  and  then  Inferring 
the  atomization  mechanlsm(s)  which  are  consistent  with  the  data.  While  not 
Inclusive,  the  proposed  atomization  mechanisms.  Figure  11,  can  be  reduced  to 
six  general  types; 

a.  Aerodynamic  Interaction  -  This  Is  similar  to  the  wind-induced  interac¬ 
tion  models  developed  for  lower  speed  jets.  Primary  breakup  Is  attributed  to 
wave  growth  on  the  liquid  surface. 

b.  Pipe  turbulence  -  Initial  jet  breakup  Is  attributed  to  Internal  fluid 
motion  produced  by  turbulence  generated  upstream  of  the  nozzle. 

c.  Wall  boundary  layer  profile  rearrangement  -  Breakup  forces  are  produced 
by  the  rearrangement  of  the  fluid  radial  velocity  profile  once  the  jet  exits 
the  nozzle  and  the  wall  boundary  condition  Is  no  longer  present.  The  redistri¬ 
bution  of  energy  results  In  radial  velocity  components  which  disrupt  the  jet. 

d.  Boundary  layer  acceleration  -  Fluid  In  the  boundary  layer  accelerates 
due  to  the  abrupt  change  In  constraint  at  the  nozzle  exit.  This  acceleration 
produces  Interface  stresses  and  surface  waves  leading  to  breakup. 

e.  Supply  pressure  oscillations  •  Jet  instability  and  cross  sectional 
variations  are  directly  related  to  variations  In  the  fluid  supply  pressure. 
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f.  Cavitation  -  Cavitation,  and  resultant  two-phase  flow,  upstream  of  the 
nozzle  exit  produce  large  pressure  disturbances  which  lead  to  atomization. 

Investigations  of  sprays  produced  by  high  pressure  liquid  jets  Injected  Into 
quiescent,  room  temperature  air  have  yielded  data  which  are  consistent  with 
the  aerodynamic  interaction  model,  if  modified  to  empirically  account  for 
nozzle  geometry  (see  Figures  12-20),  However,  the  role  of  nozzle  geometry  is 
not  well  understood  and  a  complete  model  for  the  atomization  process  could 
well  involve  components  of  cavitation  and/or  boundary  layer  effects.  It 
should  be  noted  that  the  reported  agreement  is  not  supported  by  direct  measure 
ment  of  the  droplet  size  distribution  in  the  atomization  region.  Such  measure 
ments  have  not  been  made.  Rather,  an  assumed  droplet  size  distribution,  based 
on  scaling  appropriate  to  the  modified  aerodynamic  breakup  model,  has  been 
used,  together  with  a  non-dilute  spray  code,  to  predict  the  droplet  size 
distribution  in  the  dilute  spray  region,  where  measurements  have  been  made. 
Clearly,  there  is  a  need  for  more  direct  methods  of  determining  atomization 
mechanisms  and  for  a  theory  capable  of  predicting  the  droplet  size  distribu¬ 
tion. 

Discussion  of  Dr.  Reitz's  presentation  produced  two  criticisms  of  the  results: 

1.  Theoretical  prediction  of  initial  drop  sizes  was  based  on  a  linearized 
stability  analysis  extended  to  interfacial  wave  amplitudes  well  beyond  the 
linear  regime, 

2.  Supporting  experimental  results  for  drop  sizes  are  based  on  measure¬ 
ments  of  small  numbers  of  drops  at  the  spray  perimeter,  where  flash  photo¬ 
graphs  could  be  taken.  The  interior  of  the  spray,  containing  most  of  the  mass 
flow  of  liquid,  could  not  be  measured. 

Dr,  Reitz  and  his  co-investigator.  Dr.  Bracco,  responded  that  their  studies 
represent  the  most  advanced  application  of  current  scientific  capability  to  a 
complex  physical  phenomenon  of  high  technological  relevance.  The  perceived 
shortcomings  of  this  research  should  stimulate  research  to  Improve  mathemati¬ 
cal  and  instrumental  capabilities. 

5,0  Measurement  Techniques 

An  overview  of  Spray  Measurement*",  Methods  and  Applications  was  given  by  Dr. 
William  Bachalo  of  Aerometrics,  Inc,  The  text  interspersed  with  Dr,  Bachalo's 
presentation  Figures  makes  additional  direct  connentary  unnecessary.  This 
presentation  discussed  generally  the  application  of  optical  and  physical  probe 
measurement  techniques  for  determining  drop  size  and  velocity  in  sprays.  Dr, 
Bachalo  concluded  that  the  phase  Doppler  approach,  which  has  been  the  subject 
of  his  research  activity,  offers  superior  capability  with  respect  to  the  range 
of  drop  sizes  measured  without  adjustments  to  optical  components  and  measure¬ 
ments  in  higli  droplet  number  densities  and  large  gas-phase  index  of  refraction 
gradients.  This  conclusion  was  contested  by  Drs.  Holve  and  Hess,  who  are 
Investigating  other  methodologies  for  drop  size  and  velocity  measurement. 

Dr.  Bachalo's  presentation  reflects  the  mainstream  of  research  activity  on 
Instrumentation  for  spray  characterization.  New  requirements  for  instrumen¬ 
tation  will  arise  in  relation  to  atomization  and  non-dilute  spray  research. 
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Some  demonstration  of  these  requirements  already  has  been  given  in  the  presen¬ 
tations  by  Drs.  Law  and  Reitz.  They  would  include  the  following: 

1.  Liquid  volume  fractions  sufficiently  high  to  discourage  the  use  of 
electromagnetic  radiation. 

2.  Interfacial  geometries  which  are  irregular  and  highly  complex. 

3.  Highly  transient  behavior. 

Also  related  to  these  requirements  is  the  need  to  acquire,  quantify  and 
analyze  information. 

Dr,  C.  W,  Kauffman  made  a  brief  presentation  on  neutron  radiography— a 
measurement  technique  which  has  been  demonstrated  successfully  under  high 
liquid  volume  iraction  conditions. 

An  overall  assessment  of  Dr.  Bachalo's  presentation  and  related  discussion 
is  that  measurement  of  drop  size  and  velocity  in  dilute  sprays  is  approaching 
credibility,  although  hampered  by  the  lack  of  an  accepted  standard  to  test 
accuracy.  Techniques  for  measuring  atomization  and  non-dilute  spray  behavior 
are  largely  unexplored. 

There  are  various  methods  for  incorporating  these  effects  depending  upon  the 
turbulence  model  employed.  Particularly  at  high  droplet  densities,  there  is  a 
lack  of  quantitative  data  to  guide  model  development. 

As  has  been  noted  in  the  above  discussion,  there  are  gaps  in  the  understanding 
of  non-dilute  spray  phenomena.  At  present,  these  gaps  are  being  filled  by 
approximations  based  on  assumptions  and  extrapr'ations,  largely  from  dilute 
spray  behavior.  Nonetheless,  computer  models,  incorporating  these  approxima¬ 
tions,  have  been  shown  to  be  consistent  with  overall  spray  behavior.  While 
this  may  be  partially  fortuitous,  it  suggests  that  the  overall  approach  is 
correct.  However,  better  diagnostic  techniques,  capable  of  determining  local 
gas  and  droplet  properties  in  non-dilute  sprays,  are  essential  to  provide 
detailed  verification  of  the  complex  models  now  in  use  and  under  development. 

6.0  Research  Needs 

During  the  last  two  sessions  of  the  meeting,  the  attendees  developed  a  list 
of  the  major  research  needs  in  the  general  area  of  sprays.  As  given,  and 
ammended,  these  are: 

1.  Diagnostics  for  liquid  scalar  properties. 

2.  Control  theory/dynamic  sensing  -  this  was  recognized  as  an 
application  of  spray  understanding. 

3.  Drop/turbulence  Interactions, 

4.  Interphase  transport  phenomena  -  two  phase,  including  Indirect 
drop-drop  interaction. 


5.  Critical  point  phenomena 

6.  Aerodynamic  drop  shattering  and  distortion. 

7.  Transient  drop  ignition/extinction, 

8.  Diagnostic  instrumentation  for  dense  sprays. 

9.  Highly  transient  spray  effects  -  not  covered  in  this  meeting. 

10.  Surface  wave  instability. 

11.  Cavitation  within  injectors. 

12.  Non-dilute  combustion  and  structure. 

13.  Spray-related  combustion  instability. 

14.  Identification  of  modeling  simplifications  appropriate  to  various 
regimes . 

15.  Mechanisms  of  production  of  stable  drops. 

16.  Selection  of  reference  atomizers  for  laboratory  studies. 

17.  Computational  methods  for  two-phase  flows. 

18.  Transverse  injection  -  injection  into  cross  flow. 

19.  Scaling  of  experiments. 

20.  Drop-drop  interaction  -  collision,  coalescence  and  shattering. 

21.  Surface  interactions  -  drop  with  wall. 

22.  Liquid  phase  reactions  -  monopropellant,  fuel /oxidizer. 

23.  Surface  tension. 

24.  Atomization  of  non-Newtonian  fluids. 

It  should  be  stressed  that  this  is  not  a  prioritized  listing.  No  attempt  was 
made  to  assign  priorities  to  Ine  various  research  needs.  However,  an  examina¬ 
tion  of  the  list  reveals  that  a  large  number  of  the  items  deal  with  aspects  of 
atomization.  This  was  probably  predictable  since,  as  was  noted  in  tr.e  review 
presentation,  atomization  is  the  least  understood  of  the  spray  phenomena. 

7.C  References 

One  of  the  preliminaries  to  the  meeting  was  a  request  for  relevant  background 
reference  publications  made  to  the  invitees.  The  response  to  this  request  was 
truly  overwhelming--enough  papers  were  suggested  to  require  weeks  of  careful 
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PHASE  interactions:  particle -laden  flow 


z 

o 

CO 

o: 

UJ 

CL 

CO 

o 

I- 

z 

UJ 

-I 

CD 

q: 

3 

K 


O 

o 

Z  2 

2  y 

h*  O 
O  CD 
UJ  cc 


LJ 


S  “ 

§£ 
Z  o 

<  K 
CC  O 


< 

UJ 


SI 


o  - 

Ul 

o  W 

z  ^ 


•*^2  1 

<  I-  w 

OC  UJ  X 
O  CD  K 


< 

CL 


6 

-J 

3 

O 

O 

s 

ol 

z 

Ul 

3i 

m 

z 


Ul 

e 

z 

Ul 

-J  o 

3  < 

CD  or 
o:  o 
3 

H  CL 

.  o 

O  Q 

2  >- 
O  CD 

o>- 

<  ^ 

OC  e 

H  W 
X  2“ 
Ul  Ul 


z 

UJ 


UJ 

CO 

< 

z 

a 

OC 

< 

UJ 


o 

z 


c  o 

S  UJ 
2  . 

<  o  ^ 
Ouj  I 
CO  O-  w 

^  o 

2  CO 

^  <41 

C  U-  ^ 

O  O  lo 


22 


tn 

>• 

^  < 

•  3E 
Ul  -1 

•  < 

Ui 

^  a: 
tn  < 

UJ 

y~  tn 

<  UJ 
O'  M 


oc  u 
o  o 
a.  -X 
tn  Ul 

z  > 
< 

OC  UJ 

►-  tn 
< 

UJ  X 
tn  a. 
< 

X  tn 
a.  3 
cc  o 

UJ  3 


»- 

H-  O 
en  u 
< 

u.  a 


>:  <  r 


Ul  Ul  < 
I-  -I  u 
u 


U.  IK  UJ 
Z  <  Q 

M  a. 


z 

o 


a. 

X 

3 

tn 

tn 

<c 


II.  Computation  of  two-phase  flow  no  more  complex 
THAN  FOR  SINGLE-PHASE  FLOW, 


Deterministic  Separated  Flow  (DSD  Model 


DISTRIBUTED  SOURCE  TERMS  DUE  TO  PARTICLES.  LagRANGIAN 
CALCULATION  DETERMINES  PARTICLE  TRAJECTORIES. 
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SIZE  AND  LIFETIME  OF  EDDIES. 

Computations  similar  to  DSF  models^  but  Honte  Carlo 

TECHNIQUES  USED  TO  FIND  STATISTICALLY  SIGNIFICANT  NUMBER 
OF  PARTICLE  TRAJECTORIES— CAUSING  ADDED  COMPUTATIONS. 
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APPENDIX  D 

Atomization 

Figures  by  R.  D.  Reitz 


REVIEW  OF  LIQUID  JET  ATOMIZATION 


Rolf  D.  Reitz 

Fluid  Mechanics  Department 
General  Motors  Research  Laboratories 


for  presentation  at 

Specialists  Meeting  on  Atomization  and  Nondilute  Sprays 
Sandia  National  Laboratories, 

Livermore,  CA. 

March  20-21,  1985. 


ABSTRACT 

In  the  atomization  regime  of  a  round  liquid  jet,  a  diverging  spray  is 
observed  immediately  at  the  nozzle  exit.  The  mechanism  that  controls 
atomization  has  not  yet  been  determined  even  although  several  have  been 
proposed.  An  evaluation  of  existing  theories  shows  that  eerodi'na.mic 
effects,  liquid  turbulence,  jet  velocity  profile  rearrangement  effects, 
and  liquid  supply  pressure  oscillations  each  cannot  alone  explain 
available  experimental  data.  However,  a  mechanism  that  combines 
liquid-gas  aerodiTiamic  interaction  with  nozzle  geometry  effects  is 
consistent  with  the  data.  But  the  specific  process  bj'  which  the  nozzle 
geometry  Influences  atomization  still  remains  to  be  identified. 
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OBJECTIVES 


•  Discuss  Spray  Measurement 
Requirements 

•  Review  Present  Capabilities 

•  Initiate  Discussion  On 
Current  And  Future 
Research  Goals 
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DEFICIENCIES 
IN  EXISTING  SPRAY  DATA 


Reliability 

Simultaneous  Size  And  Velocity 
Measurements 

Angle  of  Trajectory 

Drop  Concentrations 

Gas  Phase  Measurements  In  Sprays 

Polydisperse  Two-Phase  Flow 
Measurements 

Drop  Drag  Coefficient 

Turbulent  Drop  Dispersion 

Spray  Flame  Measurements 
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INTRODUCTIOII 


The  measurement  of  particle  or  drop  size  as  well  as  velocity  is 
needed  in  a  number  of  applications  associated  with  the  laser 
Doppler  velocimeter  (LDV).  It  is  well  known  that  the  particles 
used  to  scatter  light  must  be  small  enough  to  adequately  respond 
to  the  flow  velocity  fluctuations  but  large  enough  to  scatter 
sufficient  light  intensities.  A  method  for  simultaneous  measure¬ 
ment  of  the  particle  size  and  velocity  can  mitigate  the  concern 
with  the  particle  generation  by  limiting  the  velocity  measure¬ 
ments  to  the  small  particles.  In  other  cases,  fluid  velocity 
measurements  are  required  in  the  pretence  of  a  dispersed  particle 
field.  Such  two-phase  flow  measurements  may  be  associated  with 
spray  drops  in  a  turbulent  gaseous  flow,  solid  particles  in  a 
slurry,  or  bubbles  in  cavitation  studies.  Several  methods 
including  signal  amplitude  discrimination  have  been  used  to  limit 
the  fluid  velocity  measurements  to  the  small  particles.  However, 
these  methods  do  not  completely  eliminate  the  signals  from  larger 
particles  because  of  the  Gaussian  intensity  distribution  of  the 
laser.  A  method  is  required  that  can  perform  simultanous 
measurements  of  the  particle  (or  bubble)  size  over  a  large  size 
range. 

There  is  also  a  great  deal  of  interest  in  measuring  the  size  and 
velocity  distributions  of  particles  (drops,  bubbles  or  solid 
spheres) •  Current  interest  in  spray  combustion  research,  for 
example,  requires  that  these  measurements  be  made  in  complex 


turbulent,  flow  environments  and  often  in  the  presence  of 
flames.  Other  applications  demand  the  determination  of  mass  flow 
rate  which  also  requires  the  measurement  of  drop  velocity. 

In  this  presentation,  a  new  method  for  obtaining  the  aforemen¬ 
tioned  data  will  be  described.  The  method  referred  to  as  the 
Phase/Doppler  Particle  Analyzer  is  similar  to  a  standard  laser 
Doppler  velocimeter  but  has  the  added  capability  for  measuring 
particle  size.  The  presentation  will  cover  the  basic  light  scat¬ 
tering  characteristics  involved,  the  description  of  the  method, 
its  evaluation  and  examples  of  data  obtained  with  the  instrument. 
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General  Ccsomcntn  On  Particle  Sizing 


There  are  some  general  considerations  to  be  made  in  selecting  a 
particle  sizing  method.  Perhaps  the  first  is  to  assess  the 
measurement  requirements  in  terms  of  size  range,  working 
distances,  spatial  and  temporal  resolution,  and  whether  the 
actual  distribution  or  simply  distribution  parameters  are 
measured.  For  example,  the  small  angle  forward  scatter  detection 
system  measures  the  Sauter  mean  diameter  (Dj2)  based  on  the  aver¬ 
age  of  the  particles  within  the  collimated  beam  path.  Single 
particle  counters  measure  the  size  of  individual  particles  to 
produce  a  direct  measurement  of  the  size  distribution.  These 
measurements  are  made  at  points  in  the  particle  field.  Imaging 
methods  obtain  size  distributions  over  a  volume  of  the  particle 
field  defined  by  the  field  of  view  and  the  depth  of  field.  The 
size  of  the  field  of  view  as  determined  by  the  resolution 
requirements  and  the  need  to  obtain  an  adequate  statistical  dis¬ 
tribution  with  a  reasonable  number  of  recordings. 

The  type  of  sample  averaging  that  occurs  during  the  measurement 
may  also  be  important.  Modern  small  angle  forward  scatter 
diffraction  systems  and  imaging  systems  csn  produce  size  distri¬ 
butions  in  an  instant.  That  is,  unsteady  flows  may  be  measured 
to  obtain  a  time  history  of  the  spray  size  distribution.  Single 
particle  counters  obtain  size  distributions  that  are  average  in 
time.  Current  systems  obtain  measurements  at  up  to  SO, 000/sec. 
However,  in  unsteady  flows  an  ensemble  average  taken  over  several 
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cycles  may  be  required.  The  Phase/Doppler  Particle  Analyzer 
(P/DPA)  records  the  time  of  arrival  of  each  measurement  so  the 
temporal  behavior  of  the  particle  field  can  be  reconstructed. 
Time  resolved  measurements  are  useful  in  Diesel  sprays,  for 
example,  wherein  the  injection  occurs  in  a  short  period  of  time 
but  with  change  spray  density  and  size  distribution. 
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CONSIDEnATIONS  IN  TECHNIQUE  SELECTION 

MEASl'REMENT  CAPABILITY 

#  LINDTS  ON  0\XRALL  SIZE  RANGE 

e  THEORETICAL  -  e  i  RESOLUTION  LINflT  ON  SMALL  DROPS 

•  AN.\L^TICAL  DESCRIPTION  -  t.g  DIFFRACTION  THEORY  >  3;i 

•  LOSS  OF  SENSITINITY'  -  t  g  DETRACTION  BY  LARGE  DROPS  IS  CONCENTRATE 
IN  THE  FORWVVRD  DIRECTION 

•  D^'N.A.NflC  R.\NGE 

•  MEASUREMENT  R\NGE  AT  ONX  INSTRUMENT  SETTING 

•  LLNHTATIONS 

•  INSTRl-MENT  RESPONSE  FUNCTION 

•  DETECT0R/PRE.\MPLIF1ER  RESPONSE 

•  S.AMPLE  VOLUME 

•  LOSS  OF  SENSITATH' 

•  SlGN.M^TO-NOISE 

•  OBSCURATION 
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Particle  Number  Density  Considerations 


Single  particle  counters  have  been  recognized  for  their  ability 
to  obtain  size  distributions  directly  with  high  resolution  and 
can  be  coupled  with  the  LDV  to  obtain  simultaneous  measurements 
of  size  and  velocity.  However*  the  earlier  arrangements  using 
forward  scattered  light  detection  were  severely  limited  by  the 
need  to  have  only  a  single  drop  in  the  sample  volume  at  a  time. 
Off-axis  light  scatter  detection  eliminated  this  problem. 

The  modern  single  particle  counters  can  obtain  measurements  in 
sprays  with  very  high  number  densities  and  are  limited  by  similar 
conditions  to  the  ensemble  averaging  techniques.  The  primary 
limit  is  the  extinction  of  the  laser  bean.  Measurements  have 
been  made  with  the  beam  extinction  as  large  as  75%. 

Multiple  particle  or  ensemble  detection  methods  (small  angle 
forward  scatter)  are  limited  by  the  need  to  have  a  sufficient 
number  of  drops  In  the  beam  to  produce  an  adequate  signal>to> 
noise  ratio,  fingle  particle  counters  do  not  suffer  from  low 
number  densities  but  the  time  required  to  obtain  s  distribution 
may  become  excessive.  Nell-designed  instruments  provide  for  very 
large  changes  in  the  sample  volume  to  overcome  this  minor  limita¬ 
tion. 

The  number  density  or  particle  flux  is  needed  in  many  applica¬ 
tions.  Both  the  site  and  velocity  are  needed  to  accurately 
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determine  this  parameter.  In  addition,  the  sample  volume  needs 
to  be  determined.  This  will  be  discussed  in  a  later  section. 
Ensemble  methods  cannot  measure  the  particle  number  density 
directly.  Often  an  extinction  measurement  is  used  to  infer  the 
number  density. 
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•  MllTIPLE  SCATTERING 

•  ilEASlTtEMENTOFAVl 
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•  COLLEC  TION  OF  DROPS 
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Size  Distribution  Type 


In  particle  field  studies  it.  is  important  to  define  and  under* 
stand  the  differences  between  the  types  cT  distributions  that  are 
acquired.  A  spatial  distribution  is  obtained  when  a  collection 
of  drops  occupying  a  given  volume  of  space  is  sampled  instantan¬ 
eously*  as  with  a  high-speed  photographic  or  holographic  means. 

As  illustrated  in  the  adjacent  figure,  the  size  distribution  is 
sensitive  to  the  relative  number  density  of  each  size  class,  N(d) 
in  particles  per  unit  volume. 

A  temporal  distribution  is  generated  by  measuring  individual 
drops  that  pass  through  a  sampling  cross-section  during  an 
interval  of  time.  Thus,  temporal  distributions  art  generally 
produced  by  collection  techniques  and  by  optical  instruments 
which  arc  capable  of  sensing  individual  particles.  The  temporal 
distribution  depends  upon  the  particle  flux. 

The  temporal  distribution  may  be  transformed  into  the  spatial  by 
dividing  the  number  of  samples  in  each  size  class  by  the  average 
velocity  of  the  particles  in  that  size  class.  Xf  all  of  the 
drops  arc  moving  at  the  same  velocity,  the  spatial  and  temporal 
distributions  are  identical.  However,  sprays  generated  by  pres¬ 
sure  atomizers,  for  example,  will  produce  differences  in  velocity 
of  an  order  of  magnitude  between  the  smallest  to  largest  drops. 
This  emphasizes  the  need  for  simultaneous  drop  size  and  velocity 
measurements* 
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SPRAY  MEASUREMENTS 
DIAGNOSTIC  METHODS 
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•  KL\TERL\L  PROBES  -  e  g  HOT-WIRE 
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particle  Sizing  By  Light  Scatter  Detection 


There  are  several  methods  available  for  obtaining  the  particle 
site  from  the  detection  of  scattered  light.  Perhaps  the  most 
obvious  means  involves  the  measurement  of  the  amplitude  or 
intensity  of  the  scattered  light.  Particles  greater  than  the 
light  wavelength  (—3  microns)  scatter  light  in  proportion  to 
their  diameter  squared.  The  scattered  intensity  of  Xnown-sized 
particles  can  be  measured  to  calibrate  the  system  to  determine 
the  collection  efficiency  and  gain  before  measuring  site 
distributions. 

When  using  laser  beams*  the  Gaussian  intensity  distribution 
presents  an  ambiguity  due  to  the  uncertainty  in  the  Incident 
intensity  on  the  particles.  The  random  trajectories  of  the 
particles  through  the  beam  presents  a  range  of  Incident 
intensity  .i4  on  the  drops.  One  method  used  to  remove  this  uncer¬ 
tainty  is  to  accept  only  signals  from  particles  that  pass  through 
the  center  of  the  beam.  This  may  be  done  with  a  pair  of 
detectors  and  specially  designed  mashing  systems.  Unfortunately* 
the  method  requires  a  very  short  t^rXing  distance. 

Another  method  invol*«es  the  use  of  two  concentric  beams*  one  with 
a  much  amaller  diameter  than  the  other  and  having  a  different 
wavelength  or  polarisation*  This  "pointer  beam"  method  uses  the 
omaller  beam  no  indicate  when  the  particles  have  passed  through 
the  center  of  the  larger  beam.  A  deconvolution  method  has  also 


137 


been  developed  which  uses  statistics!  methods  to  resolve  the 
ambiguity* 

In  general,  scattering  intensity  methods  require  frequent  cali> 
brations,  are  subject  to  errors  resulting  from  laser  beam  atten¬ 
uation  and  improper  alignment.  High  particle  number  densities 
found  in  sprays,  for  example,  can  produce  beam  extinctions  as 
high  as  701.  Refractive  index  fluctuations  in  spray  flames  can 
produce  beam  spreading  and  deflections.  Each  of  these  effects 
produces  erroneous  measurements  resulting  in  size  measurements 
that  are  too  small. 

In  the  forward  scatter  direction,  the  angular  distribution  of  the 
scattered  light  Is  inversely  proportional  to  the  particle  size. 
The  angular  distribution  of  the  scattered  light  can  be  utilized 
to  size  particles  from  0.1  to  300  micrometers  in  diameter.  Small 
angle  forward  scattering  (ensemble  method)  is  used  for  particles 
in  the  size  range  of  S  to  300  micrometers  with  relatively  good 
accuracy.  Tor  larger  sizes  the  scattered  light  is  concentrated 
along  the  transmitted  beam  and  its  measurement  becomes  very 
sensitive  to  the  angular  resolution  of  the  detertor. 

More  recently,  the  theoretical  description  for  the  dual  beam 
light  scattering  has  been  developed.  The  light  scattering  inter¬ 
ferometry  method  bases  the  measurement  of  the  particle  size  on 
the  relative  phase  shift  of  light  scattered  from  one  beam  with 
respect  to  the  other*  This  method  hat  the  advantage  of  having  a 
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linear  reponae  to  the  particle  sizCf  ia  independent  of  the  acat 
tered  intenaity*  and  la  eaaily  combined  with  the  laaer  Doppler 
velocimeter. 
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Light  Scattering  Thieory 


Without  going  into  details,  the  light  scattered  by  homogeneous 
spherical  particles  of  arbitray  size  is  described  exactly  by  the 
Lorenz-Mie  theory.  Unfortunately,  this  method  is  cumbersome  for 
general  use  and  is  not  amenable  to  gaining  insights  into  the 
scattering  phenomena.  For  very  small  particles  (d<<\}  the 
Rayleigh  theory  may  be  used  as  a  good  approximation.  Light 
scattering  by  particles  much  larger  than  the  light  wavelength  may 
be  accurately  described  using  geometrical  optics  theory.  The 
diffraction  and  geometrical  optics  theories  are  asymptotic 
approximations  to  the  Lorcnz-Mie  theory.  These  methods  have  the 
advantage  of  making  it  easy  to  understand  the  scattering 
mechanisms. 
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CeoiBctrical  Optics  Theory 


The  geometrical  optics  description  of  light  scattering  is 
relatively  simple.  Laws  of  reflection  and  refraction  are  used  to 
describe  the  direction  and  spread  of  the  emergent  light. 

On  the  adjacent  page  a  schematic  of  the  light  scattering 
phenomena  is  shown.  A  spherical  particle  is  shown  immersed  in  a 
laser  beam.  The  lower  enlarged  figure  shows  a  representative 
pencil  of  light  incident  on  the  sphere.  Part  of  the  light  in  the 
pencil  is  reflected  and  part  is  transmitted.  The  ratio  in  the 
two  components  can  be  determined  from  the  Fresnel  reflection 
coefficients.  The  light  is  deflected  in  accord  with  Snell's  law 
and  the  laws  of  reflection. 

Using  these  basic  laws  the  scattering  coefficients  can  be 
calculated  and  are  shown  on  the  following  pages.  The  angle  6  ia 
taken  with  respect  to  the  transmitted  laser  beam.  After  obtain¬ 
ing  the  scattering  amplitude  coef ficieints  they  are  summed  and 
the  square  of  the  amplitude  is  the  intensity  (remember  that  the 
amplitude  coef ficieints  are  complex  numbers). 

Plots  of  the  scattering  intensity  obtained  from  the  simple 
theories  of  diffraction,  reflection  and  refraction  are  compared 
to  the  exact  Lorenz-Hie  theory  calculations  for  a  5  micrometer 
sphere.  The  results  Including  the  high  frequency  resonances  are 
in  very  good  agreement. 
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DR.W  INC  SHOWING  “niE  DROP  IN  THE  LASER  BEAM  AND  THE  DI\XRCENCE 
or  AN  INCIDENT  PENCIL  OF  UCHT 


iUI* 

■CAM 


T 

-L 


»A 


EMERGENT  UCHT  U>READS  INTO  AN  AKEA 


•  REFLECTION  (PERPENDICUUR  POIARIUTION) 


•  CONMENTS 

•  NOT  AFFECT  THE  ANCllAR  DISTRIBITION  EXCEPT  TKROWGH 
THE  EXPONENTIAL  PHASE  TERM 

•  SCATTERED  INTENSm'  IS  DEPENDENT  ITON  RETRACTnE  INDEX  m 

•  E.>»ONENTl  AL  PH  ASE  TERM  1$  INDEPENDENT  OF  THE  REFIUrnNE  INDEX 
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•  REFRACTION  (PERPENDICULAR  POLARIZATION) 


•  COMMENTS 

•  i  DOES  NOT  AFFECT  THE  ANGIT-AR  DISTRIBUTION  (EXCEPT  THROUGH 
THE  EXPONTNTIAL  PHASE  TERM) 

•  SCATTERED  INTENSm*  IS  DEPENDENT  UPON  REFRACTI\X  INDEX,  m 

•  EXPONXNTLAL  PRASE  TERM  IS  DEPENDENT  UPON  THE  REFRACTINX  INDEX 

•SCATTERED  LIGHT  INTENSm' 


5(0. m,#) «  Srf,/,(o.tf)4  S*>»(o,m.#)+  S‘2»(o.»n.f) 
j(o.m.#)K|5(o,in,l)  |* 
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•  EXXMPLE  or  SCATTERING  DIAGRAM 


•  COMMENTS; 


•  GEOMETRICAL  OPTICS  IS  IN  GOOD  AGREEMENT  WITH  EXACT  MIE  THEORY 
FOR  itim  DROP,  a  «  y  SS 

•  OSCE.UTIONS  IN  THE  INTENSI'H*  IS  DUE  TO  INTERFERENCE  BETWEEN 
THE  SCATTERING  MECHANISMS  (REFRACTION,  REFLECTION,  AND  DIF- 
TRACTION) 

AltOMHtlCS.lNC. 
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^  ^  -r  ^  ^  ^ 


•  COMMENTS 

•  INCIDENT  LASER  LIGHT,  Iq,  HAS  GAUSSIAN  INTENSITY  DISTRIBUTION 

•  INCIDENT  INTENSrn'  ON  DROP  DEFENDS  ON  TRAJECTORY  THROUGH 
BEAM 

•  FORWARD  SCATTERED  LOBE  -  DIFFRACTION 

•  AMPLITUDE  PROPORTIONAL  TO 

•  ANGIX.\R  DISTRIBUTION  PROPORTIONAL  TO  1/d 

•  INDEPENDENT  OF  PARTICLE  COMPOSITION 

•  -OFF-AXIS"  SCATTERED  LIGHT  -  REFLECTION  AND  REFRACTION 

•  AMPLITIDE  PROPORTIONAL  TO 

•  ANGIX.AR  DISTRIBUTION  (A\'ERAGED)  INDEPENDENT  OF  d 
m  DEPENDS  ON  REFRACTIND  INDEX 
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•  A>’GIT.\R  DISTRmUTION  OF  SCATTERED  LIGHT 

•  SMALL  ANGLE  FORWARD  SCATTER  DETECTION 

•  OPTICS 


'scat 


•  SCHENUTIC  DUGRAM  OF  ENSE\ffiLE  OF  SCATTERED  LIGHT  INTENSm- 
DISTRIBUTIONS 

•  COMMENTS 

•  AVTRAGE  OF  ALL  PARTICLES  IN  THE  BEAM  DURING  SAMPLE  ACQllSITION 

•  SENSITIVE  TO  ANGULAR  DISTRIBUTION  OF  UGHT  (PROPORTIONAL 
TO  1/rf) 

•  DECREASING  RESOLUTION  WITH  INCREASE  IN  PARTICLE  DIAMETER 

AROMETRICS.  INC. 
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COMMEyrS  ON  LIGHT  SCATTER 
DECONVOLUTION  METHOD 

•  Sensitive  To  Dcam  Attenuation 

•  Alignment  Is  Critical 

•  No  Size- Velocity  Correlations 

•  Poor  Sensitivity  To  Larger  Drops 
With  Near-Forward  Scatter 
Detection 

•  Not  Evaluated  In  Dense  Sprays 


Light  Scattering  Interferometry 


The  phase  of  the  light  transmitted  or  reflected  from  the 
spherical  particle  is  the  key  information  used  in  the  Phase/ 
Doppler  method.  Relative  phase  shifts  are  introduced  to  the 
light  wave  as  a  result  of  the  optical  path  differences  of  the 
rays  ccattered  by  the  sphere*  The  adjoining  figure  shows  a 
transmitted  ray  along  with  the  imaginary  reference  ray  deflected 
at  the  center  of  the  sphere  as  if  the  sphere  was  not  present. 

The  relationship  for  the  phase  shift  between  the  actual  and 
reference  rays  is  given.  The  t  and  t'  are  incident  and  refracted 
angles  taken  with  respect  to  the  surface  tangent,  p  is  the  ray 
of  interest  (p  >  o  first  surface  reflection*  p  ■  1  single  trans> 
mission,  p  2  internal  reflection,  etc.). 

«  PR4SE  SHIFT  DIX  TO  LENGTH  OF  OPTICAL  PATH 


A.  ■■  ^buttkih 

^  t »  fwi  lii  f*| 

4  am  drop  diamrifi 

b  m  lifhl  »t«tltO|tb 

t  V  Uits  «iib  rnprri  to  M  ima|iDary 

ro)  drS^titd  *1  tbr  tnxtt  of  tbr  tpbn* 
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The  light  scattering  interferometry  method  requires  the  same 
treansmitter  optics  as  a  dual  beam  LDV.  When  two  intersecting 
beams  are  incident  upon  a  particle,  the  particle  scatters  light 
from  each  beam  independently. 

If  a  representative  ray  from  each  beam  is  considered,  it  is  easy 
to  see  that  they  reach  a  point  P  by  different  optical  paths. 

This  is  true  because  they  enter  the  sphere  at  different  angles. 
The  optical  path  difference  results  in  a  relative  phase  shift 
between  the  rays  arriving  at  point  P.  The  relative  phase  differ¬ 
ence  can  be  computed  from  the  relationship  described  earlier.  If 
the  computations  are  carried  out  for  each  ray  incident  upon  the 
sphere,  a  scattered  interference  fringe  pattern  can  be  generated. 


•  SCKE>UTIC  or  DI  AL  SCAM  UCMTICATTDtl\C 


TWO  INTERSECTING  LASER  BEAMS  INCIDENT  UPON  THE  DROP 


OBSERVATIONS 

RAYS  FROM  BEAMS  J  AND  2  REACH  A  COMMON  POINT  F  E  <  DIFFERENT 
OPTICAL  PATHS  THAT  DEPEND  ON  THE  BEAM  INTERSECTION  ANGLE. 

AS  WELL  AS  REFRACnAX  INDEX,  m 
THE  OPTICAL  PATH  DIFFERENCE  PRODUCES  A  REUTIVX  PHASE  SHIFT 
AT/» 
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PIUSE  DIFFERENCE  BETWEEN  RAY  I  AND  RAY  2  FROM  BEAMS  1  ANT)  J 

► 

^  =  Ai  -  Aj 

^  'I  -  ^2)  “  pm(iin  r|  -  ilo  Tj)j 

PRODUCES  AN  INTERFERENCE  FRINGE  PATTERN 

COMMENTS 

PIUSE  0  IS  PROPOmTONAL  TO  THE  DROP  DUMETER, 

MEASl'REMENT  OF  ^  IMPLIES  MEASURING  THE  INTERFERENCE  FRINGE 
PATTERN 

THE  INTERFERENCE  FRINGE  PATTERN  PRODUCED  BY  A  MOMNG  DROP 
APPEARS  TO  MONT  AT  THE  DOPPLER  DIFFERENCE  FREQlTNn*. 


AE»OMnilCS,  INC. 


PHOTOGRAPH  OF  THE  SCATTERED  INTERFERENCE  FRINGE  PATTERN 


A  photograph  of  the  actual  interference  fringe  pattern  is  shown. 
The  fringe  pattern  was  produced  by  a  stationary  drop  held  at  the 
beam  intersection.  Note  that  the  fringes  are  hyperbolic  curves 
with  straight  fringes  throuah  the  bisector  of  the  beams.  The 
spacing  of  the  fringe  pattern  is  inversely  proportional  to  the 
drop  size  and  also  depends  on  where  the  observation  is  made.  The 
theoretical  predictions  provide  all  of  the  necessary  parametric 
information  for  obtaining  the  drop  site. 


AEffOVntICS.INC. 


Phase/Doppler  Particle  Analyser 


2n  order  to  obtain  drop  sire  measurements*  the  spacing  of  the 
interference  fringes  must  be  measured  accurately.  Particles 
moving  through  the  laser  beam  intersection  scatter  interference 
fringes  that  move  at  the  Doppler  difference  frequency.  The  laser 
Doppler  velocimeter  measures  this  frequency  to  obtain  the 
particle  velocity.  A  means  is  required  to  measure  the  fringe 
spacing  at  the  same  time. 

If  pairs  of  detectors  are  locates  in  the  fiold  of  the  scattered 
fringe  pattern,  each  detector  will  produce  e  Doppler  borst  signss 
but  With  a  relative  phase  shift  between  the  signals.  The  phase 
shift  between  the  Doppler  burst  signals  is  inversely  proportional 
to  the  fringe  spacing  and  hence,  directly  proportionel  to  the 
particle  site. 

Three  detectors  are  required  to  avoid  ambiguities  associated  with 
phase  angle  measurement  of  greater  than  360^.  Th\  additional 
phase  measurements  also  serve  to  extend  the  site  range  of  the 
instrument  at  one  optical  setting.  2n  addition,  comparisons  of 
the  phase  angle  measurements  between  pairs  of  detectors  are  used 
to  evaluate  the  measurements  and  to  reject  spurious  signals. 
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rH4SE/D0PrLER  STRAY  ANALYZIS 


•  OBSERVATIONS 

•  TE  ilTOR  AL  rREQlXNO  OF  THE  SCATTERED  FRINGE  PATTXRMS  UNXARLY 
RELATED  TO  PARTICLE  ATLOCITY-  #  g  LASER  DOPPLER  VTLOCIMETER 

Vflocit)  M  I  X  /d 

I  -  1  fiCAX*  INTERSECTION  ANGLE 

ft,  -  DOPPLER  DITFERCNCE  FREQITNO* 


•  SPATIAL  FREQIXNCF  OF  TH  SCATTERED  FRINGE  PATTERN  IS  LINEARLY 
REL  ATED  TO  PARTICLE  (SPHERE)  DIAMETER 


•  PRINCIPLE  OF  MEASUREAIENT  >  PAIRS  OF  DETECTORS  LOCATED  IN  THE 
FIELD  OF  THE  SCATTERED  FRINGE  PATTERN  WILL  PRODUCE  SIMILAR  DOPPLER 
BURST  SIGNALS  BUT  WITH  A  PHASE  SKRTT.  THE  PHASE  ANCLE  BETWEEN 
THE  SIGN  AL  PAIRS  IS  PROPORTIONAL  TO  THE  DETECTOR  SPACINC  AND 
THE  SCATTERED  FRINGE  SPnCINC. 


mMoiMmKM  iiinMonottAki 

•iTMOCTICTM 
AMHTUlllllMOim 


MNMTtCt.MC 


•  OPTICAL  SYSTEM 


MEASUREMENT 


DtT  1 


OCT? 


OCT  9 


f  iitciicD  ooffLCM  tunrr  sionau 


Phase  Variation  With  Drop  Sice 


The  theoretical  prediction  of  the  phase  angle  versus  the 
dimensionless  drop  size  is  shown  on  the  adjacent  figure.  The 
line  with  the  steeper  slope  identified  as  0j_3  corresponds  to  the 
detectors  with  the  greatest  spacing  while  the  lower  slope  of  ^1^2 
corresponds  to  the  smaller  detector  spacing.  In  this  particular 
example  the  ratio  of  the  detector  spacings  was  three.  With  this 
detector  arrangement*  the  signal  processor  has  sufficient 
resolution  to  provide  a  size  range  of  105  at  a  single  optical 
setting. 

The  relationships  for  all  of  the  optical  parameters  involved  were 
tested  using  monodisperssd  drops.  In  all  cases*  the  monodispcrse 
drops  could  be  measured  to  within  a  few  percent  of  the  predicted 
size. 


•  THEOR£TICAL  ANALYSIS 

«  PnOMDE  MATHEMATICAL  DESCRIPTION  OF  INTCRTERENCE  PATTERN 

•  DEFINT  PARAMETRIC  EFFECTS 

•  PRODl  CE  Fl'NCTIONAL  REUTIONSHIP  FOR  DROP  SIZE 

•  DESCRIBE  THE  SAMPU  VOLUME  AS  A  FUNCTION  OF  DROP  SIZE 


THEORETICAI.  PREDICTION  SHO\MNG  THE  PHASE  VARUTION  WITH 
DIMENSIONLESS  DROP  SIZE 


•  COMMENTS 
•  LINTAR  RELATIONSHIP 


1  -  LASER  BEAM  INTERSECTION  ANGLE 

•  THREE  OR  MORE  DETECTORS  REQITRED  TO; 

•  PRENTNT  PHASE  A^iBIGUI*^' 

•  INCREASE  DWAMIC  RANGE  AND  MAINTAIN  SENSITIYIT. 

•  REJECT  SPURIOUS  SIGN.ALS 

•  DETERNflNX  DIRECTION  OF  FRINGE  MOTIONS 


ABroMrntics.  me 
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•  \TRiriCATIONS  or  THE  THEORETICAL  PREDICTIONS 


•  SIZE  RANGE  SELECTION 

•  BEAM  INTERSECTION  ANGLE 

•  DETECTOR  SPACING 

•  OPTIC.Al  PARAMETERS  ATFECTING  MEASUREMENT 


m  <■>  rehartive  iadrx  of  drop 
X  >  laser  litht  waveleoftb 
^  •  laser  beam  ioterserlbo  ao|Ie 

#  -  observatioe  an{le 

Jt  -  obtervatioD  distaoee 
5  -  defector  tpariit| 

•  PARAMETERS  CHANCE  THE  SLOPE  (SIZE  SCALE)  OF  THE  RESPONSE  CURVX 
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Sample  Volume  Normal ieatione 


All  single  particle  counters  using  lasers  produce  distributions 
that  are  biased  due  to  the  Gaussian  beam  intensity  distribution. 
The  bias  occurs  as  a  result  of  the  larger  particles  producing  a 
detectable  signal  from  a  larger  area  of  the  sample  volume  than 
the  small  particles.  Small  particles  may  only  be  detected  when 
they  pass  through  the  central  high  intensity  region  of  the  beam. 

The  bias  may  be  removed  by  accounting  for  the  change  in  sampling 
cross-section  with  particle  site.  Analytical  descriptions  have 
been  used  by  assuming  an  ideal  Guaussian  beam  intensity  distri¬ 
bution  and  that  the  beam  intensity  is  not  affected  by  windows  or 
the  particle  field.  This  obviously  cannot  be  relied  upon  in  most 
practical  applications. 

The  Phase/Doppler  Particle  Analyser  has  incorporated  a  method  for 
measuring  the  sample  volume  cross-section  for  each  size  class. 
Since  the  sample  volume  is  measured  with  every  size  distribution, 
the  effects  of  the  measurement  environment  are  accounted  for  in 
the  sampling  statistics.  This  method  also  provides  greater  con¬ 
fidence  in  the  determination  of  the  mass  flowrate. 
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Systea  EvAluatlon 


Extensive  experimentation  was  conducted  to  verify  the  performance 
of  the  P/DPA.  The  early  experiments  involved  the  measurement  of 
monodisperse  droplets  with  the  beams  passing  through  sprays*  The 
opray  had  little  or  no  affect  on  the  measurements.  Comparisons 
of  spray  measurements  were  made  with  other  instruments.  The 
adjoining  figure  shows  an  example  of  results  obtained  by  Delavan 
Inc.  These  data  were  obtained  using  a  light  scattering  probe  for 
the  small  particles  and  an  imaging  system  for  the  larger 
particles.  The  relatively  large  dynamic  range  of  the  P/DPA 
allowed  the  coverage  of  the  entire  size  distirbution  at  one 
optical  setting* 


Other  measurements  were  made  to  characterize  the  axial  and  radial 
size  distributions  of  pressure  and  two-fluid  atomizers.  Repre¬ 
sentative  results  are  show  on  the  attached  figures.  Rote  that 
the  white  dots  within  the  histogram  bars  represent  the  data 
before  being  corrected  for  the  effects  of  the  Gaussian  beam. 


Size-velocity  correlation  data  are  also  presented.  These  data 
represent  the  mean  and  rras  velocities  for  each  particle  size 
class.  The  evolution  of  the  spray  as  a  result  of  the  relative 
velocities  of  the  particles  and  their  respective  rates  of  relaxa¬ 
tion  can  be  derived  from  these  data.  For  example*  the  size 
velocity  correlations  have  been  presented  to  show  the  spray 
development  with  axial  and  radial  distance. 
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•  PRESSURE  ATO^02ER 

•  DELAVAN  45B.  .6gpb  AT  60  p«f 

•  MODER.^TE  DROP  NT^MBER  DENSITY 


•  CONDIENTS 

•  COMPARISON  TO  DELAVA!';  INC  MEASIHEMENT 

•  AGREEMENT  ON  OiO  AND  D^i 

•  CANNOT  CONCLUDE  ACCURACY  FROM  Dj2  COMPARISON 

•  SIZE  DISTRIBUTION  GREATER  THAN  FACTOR  OF  30 

AEtOMHtlCS.INC. 
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- AtiifUMttKiCS.  INC. 


UiM 


2.3  41.1 

Pianeter  (nicroHe ters) 


’  ‘  ^  -  P  ^ 

r  -“T  ''r* s.  ’• 


Arithftetic  K?an  (P19):  15 
Area  Kean  ((20):  20 
UoluKe  Kean  ((39):  25 
Sauter  Kean  ((32):  38 

Corrected  Count;  11943 


81.7  126.0 

(ianeter  (Micrci^eters) 


flritketic  Kean  ((10):  25 
Area  Kean  ((20):  32 
VoIuKe  Kean  ((30):  38 
huter  Kean  (^2):  55 


owters)  Corrected  Count:  10335 


Drop  Size  Distribution 


Tenporal  (istrilution  — 


1,  Spatial  (istrilutien  — 


Pianettr  (si crone ters) 


liiMter  (nicroMters) 


Tempore!  and  Spatial  Size  Distributions  with 
Corrasponding  Size  -  Velocity  Correlations 
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9.0m 

I  •  10  c».  Radial  Drop  Siza  - 


Velocity  Correlations 
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lo  summary,  tbr  reroguiird  poteotial  claracteristics  of  tbe  mrtbod  m 

•  linear  relationship  between  tbe  measured  phase  an(!e  and  drop  lixe 

•  size  range  of  30  or  greater  at  a  single  optical  setting 

•  overall  size  range  of  3  to  2000  microns 

•  simultaneous  size  and  velocity  measurements 

•  relative  insensitivity  to  beam  or  light  scatter  attenuation 

•  high  spatial  resolution 

•  operition  is  similar  to  an  LDV 

•  adaptable  to  existing  LDV  systems 

•  ean  distinguish  between  gas  phase  and  droplets 

•  reduced  sensitivity  to  misalignment 

•  can  perform  measurements  independent  of  refractive  index. 

A  patent  application  has  been  Sled. 
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RECENT  DEVELOPMENT 


IN  MEASUREMENTS  TECHNIQUES 


•  Simultaneous  Measurements  of 
Size  And  Velocity 

•  Point  Measurements  In  Dense 
Sprays 

•  Minimized  Impact  of  Measurement 
Environment 

•  Improved  Dynamic  Ran^e 

•  Extended  Size  Range 
0^  to  3000  urn 


•  Spray  Flame  Measurements 


PROBLEM  AREAS 


RESEARCH  OPPORTUMTIES 

•  Verification  of  Performance 

•  Near-Nozzic  Measurements 

•  High  Number  Density  Environments 

•  Gas  Velocity  Measurements  In  Sprays 

•  Nonuniform  Seeding 

•  Spray  Flame  Measurements 

•  Drop  Data  In  Lagrangian  Frame 

•  Measurement  ol  Turbulent 
Dispersion  of  Drops 

•  Communication  With  Modellers 

•  Convincing  Funding  Agencies 
Basic  Data  Is  Unavailable 


